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Abstract 

The L2-minimax risk in Sobolev classes of densities with non-integer smoothness index 
is shown to have an analog form to that in integer Sobolev classes. To this end, the 
notion of Sobolev classes is generalized to fractional derivatives of order (3 G K + . A 
minimax kernel density estimator for such a classes is found. Although there exists no 
corresponding proof in the literature so far, the result of this article was used implicitly 
in numerous papers. A certain necessity that this gap had to be filled, can thus not be 
denied. 
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1 Introduction 

When trying to describe the goodness of an estimator, minimax performance is one optimality 
criterium possible to be consulted. The minimax risk of density estimators can be regarded 
in various settings, e.g. we differentiate between the local risk in a single point and the 
integrated risk over the whole curve. Several loss functions have been under consideration, 
such as absolute, quadratic and supremum norm, Hellinger and Kullback-Leibler distance. 
But exact asymptotics is up to now limited to a few special cases: to supremum risk in Holder 
classes, and to mean integrated square error (MISE) in analytical and in Sobolev classes. 
The latter has been examined for quite a while since in 1983, Efroimovich and Pinsker com- 
pleted the asymptotic minimax rate of the lower bound of MISE (Samarow [HJ) by the still 
lacking asymptotically exact constant, using tools that are common in information theory. 
The results were enhanced and new methods of proof found by Golubev [I] and [3], Golubev, 
Levit [§] and Schipper JOJ- Sobolev classes are classes of Z/2-integrable functions, in the 
present problem densities, for which smoothness is measured through the L2-norm of their 
/? th derivative, e N. 



s P (L) = {fei*\[ (/<*>(*)) 



dx < L J> , L < oo ( 1 ) 



Nowadays it is well known that 

2/3 



inf sup nWT E f \\f n -f\\l = 7 (A L) (l + o(l)) (2) 

fn feSfi(L) V ' 
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where 7 (ftL) = (2g + 1) ^ + M + '>)' t W 

is Pinsker's constant. Estimators attaining minimax rates of convergernce have been studied 
in abundance, e.g. kernel estimators, but also wavelet estimators and a wide range of others. 
More care has to be taken when envisaging asymptotically exact minimax estimators. 
However, the characterization of the smoothness of a given density function is incomplete 
when just assigning it to some Sp(L), £ N. Recalling the Sobolev criterion, 

L > J [fW{x)fdx = -L J |a//(a;)| 2 cL,, 

we immediately observe that Sp{L) contains densities which do not lie in Sp+i, although for 
suitably chosen V < oo and e < 1, they certainly do satisfy f \u)P +e f (u)\ 2 dw < V . 
The present article is interested in the question of whether the minimax risk can also be 
calculated for such generalized Sobolev classes. Corresponding claims are implicit in a number 
of recent papers, yet their proofs cover but the entire case. For our purpose we will employ 
the concept of the so-called fractional derivative after Riemann and Liouville, thoroughly 
discussed in Samko 

with \x\ the smallest integer greater than the positive real number x. For (3 £ N, f^> is 
the f3 derivative of /, for j3 £ R + \N it is the fi^ 1 fractional derivative of / (Samko [5], p. 
137). In case fU 3 ' is continuous and Li-integrable, then f(P\u) = (— icj) /3 /(cj). The other 
way around, if (— i • id)' 3 / is L\- or L2-integrable, the inverse transform from the Fourier into 
the time domain exists and for our purpose we define: 

fW(x) := ±J {-iufj{u)e~^dx (3) 

Existence and uniqueness of the (3^ fractional derivative of / follow thus from j- J \u)@ f(u>)\ 2 du 
< L, and Parseval's equality gives J(f^(x)) 2 dx = n= J \oj@ f (oj)\ 2 doj . 

Adopting the idea of Schipper , we find upper and lower bounds for the asymptotic 
minimax risk in Sp{L), (3 > 1/2, which are then shown to converge towards each other. 
Thereby it will be verified that the minimax risk is determined by n 2 ^^ 2 ^ +1 ^j(P,L), where 
7(/3, L) is an analogue of Pinsker's constant. A minimax kernel function for kernel density 
estimation is obtained as a byproduct from the calculation. On benefit of our a statement, it is 
for instance possible to show the asymptotically exact minimax-adaptivity of non-parametric 
estimation procedures such as the recently proposed Stein's blockwise estimator for densities 
(Rigollet [7j) and the cross-validation kernel choice for density estimation (Dalelane PP). 
The calculation of the upper bound in Schipper |10j actually holds for both entire and non- 
entire smoothness indeces, so Schipper 's Theorem 3 (Math. Meth. of Statistics (1996), Vol. 
5 No. 3, page 258-260) applies directly. To show the lower bound in Section 2, we replace 
original problem of estimating a curve by the problem of estimating a finite-dimensional 
parameter 9 (of increasing dimension) . A lower bound for the risk of such an estimator may 
be found by means of the van Trees inequality. The Bayesian risk over a least favorable 
parametric family of densities J-~q, and a least favorable prior distribution A on the space of 
finite-dimensional parameters 0, such that fe £ Sp{L) with a high probability, provides us 
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with a lower bound for the minimax risk on Sp(L). It is exactly this gap in the literature: 
fg asymptotically in Sq{L) for (3 £ N, which we have been able to close in the present 
paper. Although the demonstrations follow in general the same lines as Schipper ^HJ) the 
least favorable family of densities had to be constructed in a different way. The proof of the 
essential property (Theorem 2) applies Riemann-Liouville calculus along with approximations 
in the Fourier domain and is not similar to Schipper jlUj . 

The result for the lower bound can be considered as a special case of the theorem in Golubev 
[3], who yields lower bounds for the quadratic risk of non-parametric estimation problems in 
a variety of elliptic density classes via Local Asymptotic Normality. Unfortunately the proof 
in Golubev [2] is heavily abbreviated (the proof of a claim corresponding to our Theorem 2 
is actually omitted) and not easy to retrace. We hope that by our detailed proof, we are able 
to somehow enlighten the complicated matters. 

2 Minimax bounds 

Let Xi, . . . , X n be i.i.d. random variables with common density function / and let f n be an 
arbitrary estimator for / depending but on the sample. 

Theorem 1 (see Schipper \l(Jf Theorem 3) Let Sq(L) be the Sobolev class of those L2- 
integrable densities, which satisfy f \uj? f(u>)\ 2 du> < L for some constants (5 > and L < 00. 
Then it holds, that 

mi sup nW^ E f \\f n -f\\l < 7 (/3,L) 

fn S (L) 

The bound is maintained by a kernel estimator with the minimax kernel Kg, that is the inverse 
Fourier transform of Kq(uj) = (l — c mm (L) • with c m i n = ^ n - L7r ( 2 / 3 + 1 )(/ 3 + 1 ) ^ 

Generally speaking, the derivation of the lower bound proceeds similarly to Schipper ^U] 
(Subsection 4.1, page 262-268). However it is not the same, and so we give a little more de- 
tail. The following steps lead to the desired result, which can partly be effected analogously 
to Schipper J01> partly new proofs had to found: 

1) Construction of a least favorable parametric family of densities J-q, proof that the elements 
of J-q are contained in an e-neighborhood of the considered Sobolev class Sg(L). Both our 
center function of J-q and our perturbation functions had to be constructed in a distinct way 
to Schipper |1U| . whereas the parameter set G is the same. The proof of our Theorem 2 is 
different to that of Schipper 's corresponding lemmata (Lemma 1 through 4). 

2) Definition of a least favorable prior distribution A on the parameter set 0, proof that 
under A the elements of J-q are contained in Sq{L) itself with high probability. This time, 
Schipper's distribution A and his proof (Lemma 5, p. 266-267) are possible to transfer to our 
context. 

3) Main approximation of the lower bound via the Bayes risk over J-q with respect to A by 
means of the van Trees inequality. Again the proof of Schipper's Proposition 2 (p. 267-268) 
resembles our demonstration. 

The problem of searching a lower bound for the minimax risk over the Sobolev class Sg{L), 
can be reverted to a parametric subset of S3. Whether the minimax risk over the subclass 



3 



coincides with the minimax risk over Sp(L), obviously depends on the difficulty of the esti- 
mation problem within the subclass. We achieve our aim using the adjacent construction: 
Let us assume (3 > 1/2 and let /o be the following density from S := f] Sp 



/„(*) := ( ^- (*+W/ 2 -4 -V2<*<!/2 (4) 
I 0, otherwise 

with c a defined such that /o is a density, and for technical reason, the constant a satisfying 
J |/o(ii>) sl "yg 2 \du) = 2ir. Since J \lo@ /o(w)| 2 duo < oo for all (3 < oo, also J \uj^ fo(uj)\duj exists 

for all (3 < oo. Let gA be the indicator function on [— A-\- 1/2, A— 1/2] times the factor 2 a-i ' 
i.e. 

5a(x) := 2l L T / M+l^^] (x) - (5) 

Then /o * is a symmetric density within 5, that takes the constant value 2 A-i on + 
1,^4 — 1], and decreases smoothly towards on [A — 1,A] and [—A, —A + 1]. In order to 
constitute a sufficiently difficult estimation problem departing from this very smooth density, 
let us add some perturbation functions to fo* gA- 



1 , COS ^ 



Mx) = { l»4w4 ^<o (6) 

These perturbations will be weighted by factors 9k, where 6 = (. . . , 0-2, 0-i, 6%, O2, ■ ■ •) is 
(asymptotically) in the set: 

A (L) := \e eR^l^m <A- 2 ^ and ^Vti^y <^ 2 l\ (7) 

The set {fg\ 9 € ©^(L)} will from now on be the family of densities under consideration. 

fe(x) := fo*gA(x)(l + Y,e k Lp k (xj), (8) 

where b(9) is the normalizing constant. We cannot prove that {fg\ 9 £ &a(L)} C Sp(L), but 
instead that for all e > there exists an A £ < 00, so that for every A > A e the following 
holds: sup ee e A(L) < L + e. 

Theorem 2 Let fg and Qa(L) be defined as above. Then, as A — ► 00: 

sup ll/flll = L + o(l) 

®a(L) 

This theorem is the main assertion of our paper. Filling the gap in the hitherto existing 
literature, it enables us to go on proving the minimax bound for non-integer Sobolev classes. 
Its cumbersome and unpleasantly lengthy proof is to be found in Section 5. 
The next step leading to the lower bound requires the definition of a prior distribution A, 
which is done accordingly to ^Uj, so as to yield a parameter 9 of finite dimension: Let e > 0, 
W > and o\ > 

A(0) = [J II *°(**)> ( 9 ) 

0<\k\<W \k\>W 



4 



where 5o(.) is the Dirac function on 0, and for \k\ < W: Xk(9k) are absolutely continuous 
densities with E9 2 = o\, Q\ < G 2 a 2 (A-f.s.) for some G < oo, and the Fisher information 
d9k < (1 + e)a k 2 (with respect to the translation group {Afe(. — u)\u G E}). 
(These conditions are satisfied, for example, by independent bounded, zero mean random 
variables akCk, \k\ < W, with |^| < G, E& = 1 and the Fisher-information of the density of 
£k smaller than 1 + e.) Let us set 

e)ra(2/? + l)(/? + l)7r \2^i 
P ) 

(10) 



w 



A,L(1_ 

7T V 



4 



AA 
n 



W 



As W grows with n — > oo, the dimension of the parameter 9 will tend to infinity, allowing 
for more and more perturbation functions ip^ in the definition of fg. At the end of Section 3 
it will be shown that a\ and W of this form approximately maximize the lower bound of the 
minimax risk for the prior distribution A. 

Since A is not supported on 0^(L), we will have to show that at least the probability of 
6 G @a(L) grows with n — > oo: 

First consider that A has a bounded support, \9k\ < Gok for \k\ < W, and else 6^ = 0. With 
the above construction of cr? and W, letting A ~ Inn, condition ^2 \9k\ < A~ 2 @ +l is fulfilled 
for n sufficiently large. Lemma 1 takes care 

Lemma 1 For the prior distribution A defined above, with W and o\ as in (10), it holds 
that for n — ► oo: 

P x (0$e A (Lj) = o(n- 1 ) 

This lemma corresponds to Lemma 5 in Schipper [10 , p. 266. Its proof is exactly the same 
(p. 266-267) and we abstain from quoting it here (also see Dalelane pQ for more details). 

Theorem 3 For L < oo, j3 > 1/2 and 7(/?, L) equal to Pinsker's constant we have: 

2/3 ~ „ 

liminf inf sup E f \\f n - f\\ 2 2 > 7 (/3,L) 

n ^°° U feS (L) 

Proof Let us at first reduce the supremum of the risk by restricting the set of density 
functions. According to Theorem 2 we know that for A ~ Inn, limA^oo{fe\9 G Qa(L)} Q 
Sp(L). 

liminf inf sup E f \\f n - f\\l > liminf inf sup E f \\f n - f e \\l (11) 
n ^°° /„ Sp{L) n ^°° f n e A {L) 

For any fixed A, we find a lower bound for the supremum over G^(L) through the Bayesian 
risk with respect to A. 

inf sup E h \\ f n - M|| > inf/ E fe \\f n - f e \\ 2 dA(9) (12) 
f n e A (L) f n Jba(l) 
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In inequality (J2DJ of the proof of Theorem 2 it will be shown that l-A~ z / 2 < b(9) < 1+A~ 3 / 2 . 
Furthermore, because of orthonormality || ^fc^fcll! = Yl ®\ ^ l^fel ^ A~ 2 @ +1 , l|28() . So we 
can derive, for all 9 G G^L): 

./o - i(l ^ 

fc^o 



II/, 



0II2 



6(0) 



< max f * g A || J[_a,A] + X] 



< 



const. 



/I 



Because the set of all densities with H/H2 < 1/A) is convex, we may in (j!3|) also restrict the 
set estimators to ||/n||| < without increasing the supremum. 



inf 

/« Je A {L) 



E 



feWJn 



fe\\l dk(0) 



inf 



\\U\\1<a^ Je A (L) 

EfeWfn 



E fe \\f n -fg\\ 2 2 dA(9) 



(13) 



> _ inf 

ll/nlll— A} 

= inf 

||/n|||<A 



/0II2 dk(Q) — — P\[9 @a(L) ) (14) 



E 



feUn 



h\\l dA(9) + o(n- 1 ) 



> miE x E fg \\f n -fg\\ 2 + oirT 1 ) 

fn 

Due to H/elll < A Q 1 and < 1 it holds in (QU) that \\f n - f e \\% < 4A Q l . In (QSJ) we 

return to the complete set of estimators. 

Since fg has bounded support, i.e. [—A, A], it is equivalent, as regards the quadratic risk, 
either to estimate the function fg in the time domain or its Fourier coefficients. (fg(0) = 1 is 
known) 



(15) 



E\Efa\\fn ~ fe\\l 



KIT 

~A 



KTT 

~A 



^, ^ z ^ (r. m - % m + (fn m - t. (=) 



.4 



1 



K7T 

T 



Re/e 



.4 



K7T 

~A 



A J 



kit 



\AJ 



A 



KIT 

~A 



(16) 



The van Trees inequality (Gill, Levit |Hj) may now be applied on every single summand. For 
technical reason, the real parts are derived with respect to 0| K |, while the imaginary ones are 
derived with respect to 9i K \. 



E\Ef e 



E\Ef 9 



Re/ n 



/ KTT 



KIT 



Re/, 



Im/e 



KTT 

~A 



1 2 



> 



> 



El 



dRefg 



KTT ' 



aim/,(^f) /90_ w 



(17) 



nE A J /e (0_| re |) + J_| K | 
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where we denote If e (6 K ) = f ^^^(x)^ ^x. I K is the "Fisher information" of A K and by 
construction < (1 + e)o"~ 2 for \k\ < W and = oo for \k\ > W, respectively. Hence all sum- 
mands with |k| > W vanish from the sum. Approximations for If e (8 K ), dRefe (^f-) /<9#| K | 
und dlmfg (^) /dO_\ K \ are available from 



Lemma 2 For A - 

l + o(l) 



oo: 



2A 



00 



2^4 



3Im/ g (*f) l + o(l) 



with o(l) independent of k and The proof is postponed to Section 5. From (|15|) completed 
by (|16|) . the van Trees approximation (|17[) and Lemma 2 we thus have: 



inf£ A £ / J/ n -/ e 



inf E\E 



2A 



El* 



(t) 



K7T 

T 



K+0 V 7 V 



^) 



> 



1 ^[aRe/,(^f)/^| K 



2.4 



«+0 



2 



+ 



El 



dlmf e {!f)/dO_ lK 



2A ^ 

o<|«;|<iy 



i+°(i) V 

2VA J 



n 



E 



^+(1+^1 

1 



+ 



nE \Ifg(0-\K\) + I-\k\ 

2VJ J 

^ + (1+^1 



n 



2A{l+£) o<^<w n + 2A ^ 



(18) 



All sums obtained from W and a 2 through (|18jl. i.e. from a prior distribution A satisfying 
Lemma 1, are thus lower bounds of the minimax risk. 

What we are searching for is a bound as large as possible, we hence maximize (j!8[) subject 
to the constraint (^f) 2/3 < (1 - e)AA 2 L, such that P(9 (£ @a(L)) = o(n~ 1 ) remains 

valid. The solution to this problem is W and o\ from 1)10(1 . The maximum in (|18|) can be 
approximated as follows: 



1 



2A(1 + e)n 
1 

A(l + e)n 
1 



E 

o<|«|<w 

E 

o<«<iy 
















At 







1-I-S-V 

wJ 



W[l + o(l) 



A(l + e)n + 1 

(2/3 + 1)(^ + 1)tt 



(2/3 + 1 

2/3 



2/3 
20+ 



X 2/3+1 ^ i 

1 + e 



i 

2/3+1 



i+o(i: 



= n"WT 7 (/3,L) (l + o(l) 
Combining (JTT1) with (|T2*)l . (|15|) . (fTH|) and (|19() . we obtain the required result: 



2/3 



hm inf inf sup n^E f \\f n -fg > 7 (f3,L) 

" In S {L) 



(19) 



□ 
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3 Remaining Proofs 



Proof of Theorem 2 For fg defined in equation (JHJ), it holds that 



Te lb 



1 



< 



6(61) 
1 



{fo * 9A) {l3) + (fo *9A^2 Gk W 



09)1 



k^O 



(fo*9A) W \\ +-gL\\(fo*9Aj2 d ^k 

2 ^ ' k+0 



03) 



b{6), || (/o * <m)^||| an d II (/o * 5,4 X) 0fc¥>fc)^ll| are then considered one by one. Remember 
definition (JfJJ): fk{x) = A" 1 / 2 cos(irk/A)I(\x\ < A) for k > 0, and the same with sine for 
A; < 0. Take first the normalizing constant b{9): 

b{6) = J A J *g A {x){l + jy k ^ k {x)\dx 
f A 

= l + y^9k / fo* 9A{x)tpk{x)dx 
k^o J ~ A 

= 1 + ^2°k f fo* 9a(x)-^= cos dx + ^2 9 k f f *g A {x)^= sin ~j 
fc>0 ^ fc<0 •'~ A * 

= 1 + Y] k f fo* 9a{x)-^= cos dx + 

= »+4rE 



fc>0 



A;vrx , f A 
cos — — dx — 2 . 

A J^^A-l 



1 



/o * 9a(x) ) cos dx 



— 2 / I — ; f o * 04 (x) ) cos — — dx 



For the second term on the right-hand side we have: 



2 



cos ■ 



< 



< 



A 



dx 



\[A{2A 



k>0 



'A{A — 1 /2) 

so that for /? > 1/2 and ^4 sufficiently large, it follows that 

1_ ,4-3/2 < < l + A~ 3/2 



-2/3+1 



(20) 



(/o * 9a)^' is integrable in L2. So instead of the L2-norm of fo * qa in the time domain, by 
Parseval's equality we may as well study the L2-norm of its Fourier transform. 



(fo*9A) W) 



(fo*9A) i/3) (x) 



dx 
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= 2~ / ^/oM<mH 



1 

2^ 



2sin(v4 - 1/2)om 



(2A - l)w 



do; 



1 



< 



< 



1 

A - 1/2 _/ 2^ 
1 \ 2 1 



ujP- 1 f (uj)sm(A-l/2)oj 



dto 



A -1/2 J 2vr 

2 



1 



A - 1/2 



ll/o l 



1 1 2 



(21) 



For > 1, clearly ||/ (/3 ^l) 2 , < oo because /o lies in 5. For 1/2 < (3 < 1 we can calculate 
11/09-1) 112 _ i 



'o 



2 = 27 / l w/3 1 /o( w )| 2 ^ < ll/olll + \ 1 2/3 — 1| x , which is also less than infinity. 



The consideration of the last and most important term ||(/o * §a • ^2 Okfk)^ 111 requires a 
little knowledge about fractional derivatives. For two sufficiently regular functions / and g, 
the Leibnitz formula takes the following form: 



i=o W 



where (^) an analogue to the binomial coefficient with natural numbers: 







(3\ 



/?(/?- 2)- •• 

i\(j3-i)\ ' ' i\ (j3-i)(J3-i- 1) ■ 



As usual, (q) = 1. Now we apply this expansion to (/o * #a ■ J2 9k<Pk)^ ■ Recall the definition: 
[x\ is the integer part of a real number x, and for x positive (as in our case) \x] := [x\ + 1. 



(/o * 9A ■ E ^ 

= ||E(-)(/o*^) w (EH 



(/3"0 



i=0 
L/3J 

£ E 

i=0 

m 

£ E 

i=0 
L/3J 

£ E 

i=0 
L/3J 

= E 



(/0*5A) (i) (Ete) (/3_ lL + l^ (i) (/0*5A) W (Ete 







(/?-<) 



max 



l(/o*5A) W |||(Ete) 



G9-0 



maxM ||/ W ||i || (E^^) 



i=r/3i 

oo 



+ IIE i ) (/o*^)«(E^ 



fc^O 



+ 



i=0 



CO I. 
o lU 



i 2A 



rlKE 6 ^ 



fc^O 



E • (w) (i) (E^*) 



(0-0 



(/3-0 



+ 



|E(f)cAW>(E*») M 

i=r/3i v 7 fc^o 



(22) 



where ||/q^ ||i is of course equal to 1 for i = and finite for z = 1, . . . , \ L f3\ . When f3 G N, then 
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(■) = for all i > \(3~\ , so there is no residual. In the next step we employ: 

2 7 



(£ 



(7) 2 
2 



Wk 

k^O ~ k^O 



£<t)' 



(23) 



for all 7, proven in (|2SJ). Furthermore for 9 A (£), E |0fc| < A" 2 ^ 1 and £ ^ (if ^ < 4A 2 L 
had been determined in 0. Therefrom we can show in (|29|) that 



£' 



A:7r 
T 



2/3-/ 



< (1 + AL)A 2 - 1 for < I < 2/3 



(24) 



Hence continuing at inequality number (|22j) : 



W) 



£ 

i=0 



ll/d II ! 



i 7 2A - 1 



Y.01 

k^O 

L/3J 



A*7T 

T 



2(/3-i) 00 



+ IIE ( •) (/o*W i) (E^ 



< 



+ IIE (?) (/o**>o (i) (£<w 



(/»-*) 



(25) 



For the residual we apply Lemma 3 to our functions. It states that for functions with support 
in [-A,A\. 



E w ) • s'"-" 

i=ffl 



Setting f fo * gA and g := OkVki we proceed at inequality number (|25|): 



fo * gA • £ #fc^fc 



09) 



L/3J 



< — 

2 ~ 2A 



+ 



i=l v 7 



ll/d 


lb 


2A 




£« 




fc^O 




n/d 


)|| 
Hi 


2A 





+ o(A) (/o*w)^) f>>| 
After having derived the claim of Theorem 2, we will show in (J3(J[) that 



(26) 



II f 

||(/o*&0<™||i < ^yf , where H/^^^Hx < 00 



(LAJ)i 



(L9J)| 



(27) 
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Furthermore <J2\°k\ < A~ w+1 , but -2(3 + 1 < 0, such that can be continued as 



fo * 9 A ■ ^2 ® kLpk 

k^O 



2 < 2.1 



i=i v 7 



p(lj8J)| 



+ o(A) 



A -1/2 



1 ^/ A - 2 (3+l 



i=i v 7 



1 ||f(L/9J)|| 
2-kWJQ ! [ 



+ (A) 27r M0 , 111 o(l) 
v ; A -1/2 w 

vTfi + (i)) + o(V 1 / 2 ) +o(l) 



llf^ll < 1 



This result in connection with (|20|) and (|21() completes Theorem 2: 

O^" 1 ) H/o^lb + >/Z(l + o(l) 
\/Z + o(l) 



08) 



Still we are left to prove the intermediate assertions (|2H|) . (|24j) and (J23 
As an exception to the ordinary case, sine and cosine enjoy an easy to calculate fractional 
derivative: sin^ 7 )(ax) = a 7 sin (ax + jir/2) and the like for cosine (Samko Hj, p. 174). 
Obviously, the orthogonality between our functions <p k is preserved through derivation. 



/(E^i 7) (*)) 2 ^ 

J k^O 



£*2 ^\x?dx 

k^O 
A 

A 



E' 

fc>0 



1 / kir\ 2 ' y 2 f ^ t7TX 1 7T \ 



A \A J 
~A) A 



cos 



— / cos 

A 



A 

kirx 77r 



7 fc<0 v 7 



1 / kir\ 2 ^ . 2 {kirx | 77r 



sin 



V A + 2 



ri.r 



/fcvr\ 27 1_ 



I/ 03 — &+ E"Ht) sX/ 1 



E< 

fc>0 

E< 

fc>0 

E^(t)^-E<^) 27 ^ 

fc>0 1 ' k<0 V 7 

E' 



A 

sm" 

A 

kirx 
~A 



2 / kirx "fir 



A 



— I dx 



tlx 



(28) 
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Referring to step pij). < Z < 2/3: 
kir\ 



E <^T + E «(f j 



2/3 / -J 

V a J 



2 f kir^ f kir\ 



< av-* y: d i+ A ~ i e 9 i{ k i) 

C\-L \ A1 l~ I 7„l^ A 2 \ / 
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0^|fc|<A 2 /7T 



|fc|>yl 2 /7I 



< A 2P-lJ2\0 k \+A- l J2 

k^O k^O 

< A 2/3 ~ l ■ A~ 2 ? +1 + A~ l ■ ALA 2 
= A 1 ^ 1 + ALA 2 ~ l 

< {l + AL)A 2 ~ l 



(t) 
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Proof of (EZJ: 



||(/o*^)^D| 



« w /o(*(w; 



2sin(A - 1/2)cj 



(2A - l)w 



do; 



1 



A - 1/2 



u L/3J /oM sin(i4-l/2)w 



< 



II/, 



A - 1/2 



exists, because we chose /o € S. This concludes the proof of Theorem 2. 



(29) 



(30) 



□ 



Lemma 3 For functions f and g, which are both L2-integrable, sufficiently regular and 
have support in [—A, A], it holds that 



E 



o(A 2 ) 



\g\\l 



Proof This proof takes a detour via Fourier coefficients. Begin with the following discus- 
sion: The power function is an analytical function. We may thus for instance expand (J^f 
into an infinite Taylor series at point iOz^hL m 



kit 
~A 



E 

i=0 



(3\ /Xiry /(k - X^xP-i 



i \AJ 



A 



We cut the Tailor expansion of 



after \_f3\ and bound the residual. 



IE^t)'^""' 

oo 

IE 

i=0 



/3 \ /Xir\W]+i /(k- X)ir\P-W]-i 



\P\+iJ \AJ 



A 
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p...(/3-\p\-i + l) / Attn / («- A)7 

1^ ([/?!+ <)! vaJ V A 

| /?•••(/?- +1) ,\^\(f]™ (i3-\P\)...(P-\P\-i + l) /Att 



< 



Avr I T/31 r ' 
A I 



+*)••• a/3i + i) 

g I (/3- [/?])•••(/?- \0\ -i + l) ^Xiry ^K-X)^ 
i=o A A 



A 



The product (ft — [/?]) ■■■(/? — [~/3] —i + l) consists of i factors, which are all negative. We 
can write \{P - ■ ■ ■ (0 - \P] - i + 1)| = (-!)'(/? - ■ ■ ■ (/3 - \p] - i + 1), such that 



1 E (?) (f )* (^f| 

//3\ AttiW^ r/?]) ■■■(/?- + i) 



Att 


i 


(k — X)tt 




A 




A 





(3 



E 

i=0 



-)' 



/ \ AvriCT / 



\n I m ( {\k-\\ - |a|)tt ^ /3- r/?! 



A 



Since we know that -1<P~\P] < 0, we can approximte ( (I^H^ )/?-^] = O(^-^l) = 
o(^4). Now we expand the tail of our Leibnitz formula into a Fourier series and plug in the 
bound of the Taylor series: 



(0 . g W-i) 



II t (!) /" 

- ^e(e(^e^)^(^)' 

- ^e(e(^e(^)*/(^)(^)"X^)' 

- ^E (^E E (!) (t)*(^)"/(t) K^j 
s ^E(^E^)i 

- °(^E^E 



/9 



A^ 



(ft — A)7T 
I 

(K — A)7T 
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o(A 2 ) ' 1 



-y 

2A ^ 

2 





fiTT 




/(f) 




(r£i) 


A 







A 



1 -^ 7r 



Aez 



2 



(K — A)7T 

A 



1 ^ 
2A^\~A 

(k — /i)7T 



A 



2A^\\A 



A 



-y\ 



2A^\\AJ J \~A 



= o{A 



i 

2A 

2 



El? 



(k — \)ir 
A 



2 / 1 v^U/ ( K ~ / i ) 7rN 2 



2A 



For growing A, the Fourier expansion approaches the Fourier transform, and hence 

2 
2 



1 1 (!) /( 







= o{A 2 
= o(A 2 ) fW 



\(3}) Utt 

2 



oM 1 + 



o(l)) 



i 



(31) 



which is the statement of Lemma 3. 



□ 



Proof of Lemma 2 We start with 



A {df e (x)/dB K f 
a fe(x) 



dx 



1 



1 



b 2 (0) J-a fe(x 

,-A 



f A I 2 

-fe(x) fo * 9A(y)<fK(y) dy + f * 9A(x)(p K (x) dx 

J-A 



= W(6) I A (/ J Q * 9A(y)<p K (y) dy^j -2 j J Q * 9a{v)<P k (v) dy f * g A (x)(p K (x) 



+ 



1 / \ i 

TJX (/o *5A(x)^ K (x)J 



dx 



b 2 (6) . 



fo * 9A(y)^K(y) dy 



b 2 {B) lL A f ° 



9A(y)<p K (y) dy 



+ 



+ 



(/o*0A(a:)) 



62( ^ } 7 - A E(k/o * <m(x) (l + E *a¥>a(z)) 

W V A^O 7 



dx 



fa*9A{x)^ 2 K {x) 
b(9) J-a 1 + £ 

A^O 
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1 



+ 



b 2 {9) 
1 

W) 



2A-1 



dx-2 



1 



l \2A - 1 



fo * 9a{x) ) <p K {x) dx 



2 A 



-a i + E 0\<p\fc 

A^O 



dx - 2 



i 



2A-1 



IA-1 



fo*9A(x)j ip 2 K {x) 
1 + £ 8\Mx) 

A^O 



dx 



The leading term is ^^vlw ^ Due to 1 ^ ^a(*)| < A' 1 / 2 £ |0 A | < ^-1/2.^-2/3+1 
< A" 1 / 2 , we know it lies in the interval ((2A - l) -1 ^ -M~ 1/2 ) _1 , (2A - 1) _1 (1 - ^ 1 / 2 )" 1 ). 
Moreover from (1201) we have 1 - ^" 3 / 2 < b(9) < 1 + A~ 3 / 2 . For A — ► oo we obtain: 



b 2 {9) 



+ O 



(A-1/2)VA 
l + o(l))(9(A- 3 ) + (l + o(l) 



+ 



1 



l + o(l) 



l + o(l) / l + o(l) \ 
2A-1 V(^-V 2 )^/ 



2A 



Y + ( A ~ 2 ) 



2A 



Refe (^) can be expressed as A 1 ! 2 J fg(x)(f\ K \(x)dx, yielding 



do 



d 



09 



A \ fe(x)<P\K\(x)dx 



b{9) 



I -A 



fe(x) fo * 9A(y)<P\ K \(y)dy <p\ K \(x)dx + / / * 5j4 (x)^ 2 K |(x)(ix 



/-A 



-A 



6(0) 



A 



fo * gA(y)v\ K \(y)dy + / /o*^^)^^^)^ 



74 



Re/(9 ( — ) " /"^ 1 f A ( 1 \ 

" y A 2A~^T v\K\(y) d y + 2 j A X 2A-1 ~f°* 9A ( y n ^\^y) d v 

j_ A 2A- 1 y H ~ 2 Sa-xK ^A- 9A ^) ^U^ dx 



m 

l + o(l) 

l + o(l) 
2V3 



O + O 



(A -1/2)VAJ J 1 + 0(1) 



+ 



/I 



+ 



2A-1 V(^-l/ 2 )^ 



A similar result is obtained for Im/e (^f ), whereby f fo * 9A(y) i P-\K,\(y)dy = simplifies the 
task, because the sine function is anti-symmetric. □ 
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